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Abstract 

In this paper, we propose a fast search-free method for direction-of-arrival (DOA) estimation with 
coprime arrays. A coprime array consists of two uniform Unear subarrays with inter-element spacings 
MX/2 and NX/ 2, where M and N are coprime integers and A is the wavelength of the signal. Because 
uniform linear arrays enjoy many processing advantages over arbitrary geometry arrays, our strategy is 
to exploit the uniform linear structure of the subarrays in coprime arrays. We first estimate the DOA's 
for each uniform linear subarray separately and then combine the two estimates from the two subarrays. 
For combining the two estimates, we propose a method that projects the estimated point in the two- 
dimensional plane onto one-dimensional line segments that correspond to the entire angular domain. By 
doing so, we avoid the search step and consequently, we greatly reduce the computational complexity 
of the method. 

Index Terms 

DOA estimation, coprime arrays, coprime sampling, uniform linear array 

I. Introduction 

Direction-of-arrival (DOA) estimation using sensor arrays is a problem that is frequently encountered 
in many engineering areas including radar, sonar, and wireless communication, and it has been studied 
for several decades. For arbitrary geometry arrays, popular DOA estimation methods include beamspan 
ni, Capon 111, MUSIC and several versions of maximum likelihood estimators |l4l|. For finding the 
optimal estimate, they all require a computationally expensive search step in a nonconvex space. For 
uniform linear arrays (ULAs), the array response vector has a Vandermonde form, which allows the 
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search step to be replaced by polynomial rooting. These so called search-free algorithms include IQML 
in, MODE (Method of Direction Estimation) Ig), root-MUSIC Q and ESPIRIT [8|. Besides, there are 
also approaches that transform arbitrary arrays to equivalent ULAs so that one can take the advantage of 
efficient algorithms for ULAs, e.g., array interpolation ||9l, manifold separation (TO'l and Fourier-Domain 
(FD) root-MUSIC 1 1 1 1 . These approaches basically try to approximate the steering vectors by virtual 
ULAs. Their drawback is that for achieving satisfactory performance, they require a large number of 
virtual arrays, which entails increased complexity. Furthermore, they suffer fixed error at high signal-to- 
noise ratios (SNRs) and thus, do not converge to the Cramer-Rao bound as the SNR increases. 

Coprime signal processing has recently been studied in |[T2l . |[T3l . |[T4ll . It has been shown that coprime 
processing allows us to sample a signal sparsely and estimate the parameters of the signal at a higher 
resolution |[T3l . The parameters can be the frequencies of temporal signals or the DOA's of spatial 
signals. A coprime array consists of multiple ULAs with different inter-element spacings that are mutually 
coprime. Although the DOA estimation with ULAs has been thoroughly studied, there has not been 
much work for designing special algorithms that take advantage of the partial uniform linear structure of 
coprime arrays. In 1 12 1, the authors apply the MUSIC algorithm to coprime processing, where they mainly 
consider the problem from the perspective of degree of freedom, i.e., how many signals the coprime array 
can identify. To achieve additional degree of freedom with coprime arrays, that is, to enhance the rank 
of the covariance matrix of the received signal, they propose the use of spatial smoothing. Nevertheless, 
the proposed method is not search-free, and its accuracy and computational complexity have not been 
investigated in detail. 

In this paper, we address the problem of DOA estimation with coprime arrays with the emphasis on 
estimation accuracy and computational complexity. We propose a search-free method that exploits the 
uniform linear structure of the subarrays of the coprime arrays. Considering that the DOA estimation 
using ULAs is fast and accurate, we first estimate the DOA for each subarray separately. Since the inter- 
element distance for each subarray is larger than half-wavelength of the signal, ambiguity is present and 
as a result, the DOA's cannot be uniquely determined. To resolve the ambiguity, we use a projection-like 
method that combines the estimates from different subarrays, where the correctness of the estimate is 
guaranteed as a consequence of the coprimeness. 

The paper is organized as follows. In Section 2, we introduce the model and briefly review the DOA 
estimation with ULAs. The proposed algorithm is discussed in detail in Section 3 and the numerical 
experiments are presented in Section 4. Section 5 provides final remarks. 

We use the following notation: (•)^ and (•)^ denote the transpose and the conjugate transpose. 
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respectively; C refers to the complex space; Sk,i signifies the Kronecker delta; E(-) stands for the 
expectation operator; and tr{A} is the trace of the matrix A. 

II. DOA ESTIMATION FOR ULA 

The problem of finding the DOAs of D narrowband plane waves impinging on a ULA of L sensors 
can be modeled as follows [15| : 

y(A:) = Ax(fc)+w(fe), k = l,---,K, (1) 

where y{k) e C^^^ is the signal received by the sensors at the kth time slot; K is the number of snapshots; 
x(A;) G C^^^ contains the complex envelopes of the emitter signals; ■w{k) is the noise process; and 

A= [ai,--- ,a^] gC^x^ (2) 

with 

a.d = [l,exp(jpiV'd), • • • ,exp{jpL-iipd)V ^ C^""^ 

d=l,---,D; (3) 

PiX/2 is the distance from the Zth sensor to the reference point; tpd is the DOA of our interest. The 
number of sources, D, is assumed to be known. Thus, ^' = [tpi, - ■ ■ ,1^0] is the vector of DOAs we 
wish to estimate. We assume the signal x(fc) and the noise w(A;) are independent zero-mean complex 
Gaussian random processes with the following moments: 

E[x(A:)x^(0] = PSk,i, E[x(fc)x'^(Z)] = 0,k^l (4) 

E[w{k)w"{l)] = a^I6k,i, E[w(fc)w^(0] = 0,k^l (5) 

with I being the identity matrix, P, the signal covariance matrix, and cr^, the noise power. The covariance 
matrix of the received signal y(A;) can be written as 

R = APA" + a^I. (6) 

We denote by R the sample covariance matrix 

K-l 
k=0 

The eigendecomposition of R in Q can be written as 

R = E,A,E« + E,A,E^, (8) 
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where € C^^-° and A„ E ,^{l-d)x{l-d) 

are diagonal matrices that contain the eigenvalues of the 
signal and the noise subspaces, respectively, whereas G C^^^ and E„ G (£Lx(L-D) composed of 
the eigenvectors of the signal and the noise subspaces, respectively. 
For ULAs, the steering vector becomes 

[1, exp(jV'i), • • • , exp(j(L - 



(9) 



There exist a large number of fast algorithms for estimating DOA's with ULAs, among which MODE is 
the leading candidate llT6l . ||T7 il. Since we later use MODE as part of our algorithm, we review it here 
briefly. The MODE algorithm was originally proposed in |6]. It estimates the DOA's via the minimization 
of the following function: 



where 



b = argmin tr{nBE^WE^} 



b = [5o,6i, • • • ,bD]^, 
Hb = B(B"B)^B^, 



1 



a 



L-D 

and BisaLx(L — D) Toeplitz matrix defined by 

ho ■■■ ho 



tr{AAr}, 



(10) 

(11) 
(12) 
(13) 
(14) 



B^ 



ho 





^0 



(15) 



••• ho ■■■ ho 

Let z = exp{jipi), and denote by b{z) the complex-valued polynomial 



b{z) = boz^ + biz 



ho 



D 



boUi 



Z - Zi 



(16) 
(17) 



The angles of the roots of the polynomial are the estimates of the DOA's. It was shown in |l6l that the 



MODE algorithm is an asymptotically efficient estimator of the DOA's. To minimize ( [T0| ), both iterative 
and noniterative steps have been proposed. See lfT6l . ifTTll for detail. 
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Fig. 1. a coprime array with two ULAs, wliere M — 7 and N — 5. 

III. DOA ESTIMATION FOR COPRIME ARRAYS 

We consider a coprime array with two uniform linear subarrays. We assume that the two subarrays 
have inter-element spacing MX/2 and NX/2, respectively with M and N being coprime. Fig. [T] shows 
the case for = 7 and M = 5. Because the two subarrays share the first sensor, the total number of 
sensors L is equal to M + — 1. The corresponding steering vector in Q becomes 



1 pj^^M j^a2M ... JMN-l)M 



(18) 



For such arrays, the efficient methods for ULAs are not directly applicable. We can of course use the 
algorithms for arbitrary arrays, but they are slow and the uniform linear structure of the subarrays is 
wasted. Our objective is to develop an algorithm that is applicable to coprime arrays while at the same 
time enjoys the efficiency brought by the uniform linear structure of the subarrays. The proposed approach 
is composed of two steps. We first estimate the DOA's from each subarray separately and then combine 
the two estimates in an innovative way. 
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Fig. 2. Ambiguity on tlie angular domain, case I. 
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Fig. 3. Ambiguity on the angular domain, case 2. 
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A. DOA estimation using a single subarray 

For a subarray, the inter-element spacing is NX/2 with > 1. This has the same effect as we 
undersample a signal by a factor of N in the temporal domain. According to the basic sampling theorem, 
we will have ambiguities, or aliasing, in the angular domain. Specifically, if we use the same algorithm 



as described above, the polynomial in ( [T7| ) becomes 

D 

b'{z) = b',ll{z''-zi). (19) 

i=l 

This polynomial has N x D roots instead of D. If Zi is a root of ( [TT] ), it is still be a root of ( 19i. The 
problem is that aliasing occurs with the period of 2it/N; Zjexp(jn^) for n = 1, • • • , — 1 are also 
the roots of ^T9\ , as shown in Fig. |2] Therefore, we have no way of distinguishing the Zi we want from 
the remaining — 1 roots. 

Note that in this step, it does not matter which algorithm is used. In the experiment, we use MODE, 
which was shown to have the best performance for ULAs ifTTl . 



B. Combination of the estimates from subarray s 

The present ambiguities in the estimates from the two subarrays notwithstanding, we show that we are 
able to disambiguate them. For example, in the case of one source and no noise, in Fig. |2j we can see 
that the estimates from the two subarrays coincide at the red arrows. Therefore, we can be very confident 
that the position of the red arrow is the DOA. In the presence of noise, however, coincidences as the one 
in Fig. [3] are unlikely, and consequently, it becomes difficult to tell the exact DOA. Furthermore, in the 
case with multiple sources, obtaining the DOA's in this way becomes impossible. 
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We propose a projection method that can uniquely determine the exact DOA. Recall that the aliasing 
period in the angular domain for the subarray is 2tt/N. Thus, we only need to take the values between 
[— vr, — 7r+27r/A^]. The ranges of the outputs of the two subarrays define a rectangle in the two-dimensional 
plane as shown in Fig. |4j where N = 2 and M = 3. It is not difficult to see that the 45 degree oblique 
line segments colored in blue correspond to the entire angular domain. The Chinese remainder theorem 
guarantees that the map is one-to-one and onto as a result of the coprimeness. This is the key point 
of our algorithm. In Fig. |4j the entire angular domain consist of four line segments; LI corresponds 
to [-7r,-7r/3]; L2, [-7r/3,0]; L3, [0,7r/3]; and L4, [7r/3,7r]. For general M and N, there will be 
M + N — 1 oblique line segments. Suppose the outputs of the first and the second subarray are ip^^'^ 
and respectively. The two estimates specify a point in the plane. We project the point onto the 
oblique line segments that corresponds to the entire angular domain. Specifically, we seek the point on 
the oblique line segments that is nearest to that point specified by the two estimates. This ensures that 
the combination is optimal and the result is valid. For point B in Fig. |4j we can simply draw a line 
that is vertical to the oblique line and measure the distance from B to the intersections. For point A, 
however, one intersection falls outside the rectangle area and should be wrapped around. As a result, 
the intersection is actually on L3 in the figure. After the nearest point on the line segments is found, it 
is easy to calculate the corresponding DOA by solving a set of modular linear equations, or it can be 
precalculated and made into lookup table. 

C. Combination of the estimates in the case of muhiple sources 

In the case of multiple sources, the method gets more complicated. Suppose that the outputs of the 
first subarray are V'l^'' and ip^'^ and the outputs of the second subarray, V'l^^ and '4)^\ We have no way 
of knowing which value corresponds to the first and which to the second target. Therefore, we have to 
check two possibilities. In Fig. |4j '4)^i \ '4'^\ V'l^^ and ^2^^ define four points on the two-dimensional 
plane, which are marked A, B, C, and D. Accordingly, we can choose for projection the estimates given 
by the pair {A,B) or {C,D). To decide which pair is more likely, we can simply choose the one with 
the minimum sum of distances to their nearest oblique lines. For the case of D targets the number of 
possible combinations is equal to the factorial of D. Although this may be very large for large D, we 
can use greedy algorithm to find the optimal combination with linear time complexity. 
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IV. Simulation 

In this section, we use a numerical experiment to justify our proposed method. We show the per- 
formance comparison between the proposed method, the FD root-MUSIC method, and the Cramer-Rao 
bound. The FD root-MUSIC uses Fourier series coefficients to approximate the null-spectrum function 
ifTTI . It was shown to have better performance than other search-free methods lITSl . 

In the simulations, we set = 7, Af = 5. For the FD root-MUSIC, 99 Fourier series coefficients 
are used for approximation. We only test the case of one source with ipi = O.lvr. Fig. |5] shows the 
mean square error (MSE) performance for different values of SNR. We can see that at low SNRs, the FD 
root-MUSIC outperforms the proposed method. At high SNRs, however, the performance of the proposed 
method approaches the CRB while the FD root-MUSIC suffers fixed error due to its approximation. In 
Fig. [6} we plotted the MSE performance for different number of snapshots. Similarly, the FD root-MUSIC 
is better when the number of snapshots is small. On the other hand, the proposed method approaches 
the CRB when K > 600 while the FD root-MUSIC does not improve when K > 300. 
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Fig. 6. Mean square error performance for different number of snapsfiots. (SNR=— 12dB) 



V. Concluding remarks 

In this paper, a fast search-free DOA estimation algorithm for coprime arrays is proposed. It exploits 
the uniform Unear structure of the subarrays by first estimating the DOA separately within each subarray 
and then combining the estimates from different subarrays. Simulation shows that the performance of the 
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proposed algorithm is close to the FD root-MUSIC method at low SNRs and much better at high SNRs. 
We point out that the proposed method has much lower complexity. 

In the combination step, we seek the point with the shortest distance under an implicit assumption that 
the error variances along the two dimensions are equal. In practice, however, this may not be the case. 
To take the different error variances into consideration, one needs to construct an oblique projection to 
find the point on the line segments. 

Here, we only discussed the case with two subarrays. For a coprime array with m> 2 subarrays, the 
state space will be m dimensional hyperrectangle, but the entire angular domain can still be mapped to 
fine segments and the same idea for estimation can be reproduced. Therefore, the presented results can 
be extended to coprime arrays with multiple subarrays. 

References 

[1] J. Mayhan and L. Niro, "Spatial spectral estimation using multiple beam antennas," Antennas and Propagation, IEEE 

Transactions on, vol. 35, no. 8, pp. 897-906, 1987. 
[2] J. Capon, "High-resolution frequency-wavenumber spectrum analysis," Proceedings of the IEEE, vol. 57, no. 8, pp. 

1408-1418, 1969. 

[3] R. Schmidt, "Multiple emitter location and signal parameter estimation," Antennas and Propagation, IEEE Transactions 

on, vol. 34, no. 3, pp. 276-280, 1986. 
[4] R Stoica and A. Nehorai, "Performance study of conditional and unconditional direction-of-arrival estimation," Acoustics, 

Speech and Signal Processing, IEEE Transactions on, vol. 38, no. 10, pp. 1783-1795, 1990. 
[5] Y. Bresler and A. Macovski, "Exact maximum likelihood parameter estimation of superimposed exponential signals in 

noise," Acoustics, Speech and Signal Processing, IEEE Transactions on, vol. 34, no. 5, pp. 1081-1089, 1986. 
[6] P. Stoica and K. C. Sharman, "Maximum likelihood methods for direction-of-arrival estimation," Acoustics, Speech and 

Signal Processing, IEEE Transactions on, vol. 38, no. 7, pp. 1132-1143, 1990. 
[7] B.D. Rao and K. V. S. Hari, "Performance analysis of root-MUSIC," Acoustics, Speech and Signal Processing, IEEE 

Transactions on, vol. 37, no. 12, pp. 1939-1949, 1989. 
[8] R. Roy and T. Kailath, "ESPRlT-estimation of signal parameters via rotational invariance techniques," Acoustics, Speech 

and Signal Processing, IEEE Transactions on, vol. 37, no. 7, pp. 984-995, 1989. 
[9] B. Friedlander and A. J. Weiss, "Direction finding using spatial smoothing with interpolated arrays," Aerospace and 

Electronic Systems, IEEE Transactions on, vol. 28, no. 2, pp. 574-587, 1992. 
[10] F. BeUoni, A. Richter, and V. Koivunen, "DoA estimation via manifold separation for arbitrary array structures," Signal 

Processing, IEEE Transactions on, vol. 55, no. 10, pp. 4800^810, 2007. 
[11] M. Rubsamen and A. B. Gershman, "Direction-of-arrival estimation for nonuniform sensor arrays: from manifold separation 

to Fourier domain MUSIC methods," Signal Processing, IEEE Transactions on, vol. 57, no. 2, pp. 588-599, 2009. 
[12] P. Pal and P. P. Vaidyanathan, "Coprime sampling and the MUSIC algorithm," in Digital Signal Processing Workshop and 

IEEE Signal Processing Education Workshop (DSP/SPE), 2011 IEEE. IEEE, 2011, pp. 289-294. 



January 18, 2013 



DRAFT 



11 



[13] P. P. Vaidyanathan and P. Pal, "Sparse sensing with co-prime samplers and arrays," Signal Processing, IEEE Transactions 

on, vol. 59, no. 2, pp. 573-586, 2011. 
[14] P. P. Vaidyanathan and P. Pal, "Theory of sparse coprime sensing in multiple dimensions," Signal Processing, IEEE 

Transactions on, vol. 59, no. 8, pp. 3592-3608, 2011. 
[15] H. Krim and M. Viberg, "Two decades of array signal processing research: the parametric approach," Signal Processing 

Magazine, IEEE, vol. 13, no. 4, pp. 67-94, 1996. 
[16] J. Li, P. Stoica, and Z. S. Liu, "Comparative study of IQML and MODE direction-of-arrival estimators," Signal Processing, 

IEEE Transactions on, vol. 46, no. 1, pp. 149-160, 1998. 
[17] H. L. Van Trees, Optimum array processing, Wiley-Interscience, 2002. 

[18] A. B. Gershman, M. Riibsamen, and M. Pesavento, "One- and two-dimensional direction-of-arrival estimation: An overview 
of search-free techniques," Signal Processing, vol. 90, no. 5, pp. 1338-1349, 2010. 



January 18, 2013 



DRAFT 



